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We propose a connection between conformal field theory (CFT) and the exact solution and in- 
tegrability of the reduced BCS model of superconductivity. The relevant CFT is given by the 
5'(/(2)a ; -WZW model in the singular limit when the level k goes to —2. This theory has to be per- 
turbed by an operator proportional to the inverse of the BCS coupling constant. Using the free field 
realization of this perturbed Wess-Zumino-Witten model, we derive the exact Richardson's wave 
function and the integrals of motion of the reduced BCS model in the saddle point approximation. 
The construction is reminiscent of the CFT approach to the Fractional Quantum Hall effect. 

PACS number: 11.25.Hf, 74.20.Z, 04.20.Jb 
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I) INTRODUCTION 

The BCS theory has been used for decades to de- 
scribe the superconducting properties of "low T c " metal- 
lic materialso. The starting point of the theory is a 
Hamiltonian which describes the attractive interaction 
between the electrons in well defined energy levels. The 
grand canonical BCS wave function gives a very accu- 
rate solution of the BCS Hamiltonian in the limit where 
the number of electrons, N e , is very large. The BCS 
wave function may also be projected to a fixed number of 
electrons, N e , giving essentially the same physics when 
N e >> 1. However for small values of N e , one has to 
use exact analytical or numerical methods to obtain re- 
liable results. The study of small fixed- iV e superconduc- 
tivity has a long story which goes back to an old ques- 
tion posed by Anderson as to what is the smaller size of 
a metallic particle to remain superconductinga. Recent 
experimentsa involving aluminium grains with nanome- 
ter size have inspired a number of theppejtical works where 
Anderson's question is reconsideredoT3. 

In this paper we shall be concerned with the exact ana- 
lytic solution of the reduced BCS Hamiltonian proposed 
by Richardson in a series of papers between 1963 and 
1977li3cJ. This solution emerged in the framework of 
Nuclear Physics and has passed unnoticed by most of the 
physics comunity until the recent upheaval in ultrasmall 
metallic grains. A closely related work is that of Cam- 
biaggio, Rivas and Saraceno (CRS) who proved recently 
the integrability of the reduced BCS Hamiltonian with- 
out recourse to Richardson's solution^. These authors 
found a set of integrals of motion whose number equals 
that of the degrees of freedom of the system. Our aim 
is to show that Richardson's solution, together with its 
integrability, can be naturally understood in the frame- 
work of conformal field theory (CFT). In other words, 
we shall propose a correspondence between CFT and the 
BCS theory which gives a neat picture of the Richard- 
son's wave function and the conserved quantities found 



by CRS, which are in that sense unified in a common 
framework. Our work is reminiscent r-tp the application 
of CFT to the Laughlin wave f-upctionEZI of the Fractional 
Quantum Hall effect (FQHE)E3. 

The organization of the paper is as follows. In sections 
II and III we present brief reviews on the exact solu- 
tion of the reduced BCS Hamiltonian and the free field 
realization of the SU (2)k-WZW model. In section IV 
we derive the exact Richardson's solution and we show 
integrability of the BCS model using the free field re- 
alization of the SU (2) k-WZW model. In section V we 
explain the BCS/CFT connection in a second quantized 
language. In section VI we explore the analogies and dif- 
ferences between the CFT approaches to BCS and the 
FQHE. Finally, in section VII we state our conclusions 
and prospects of future work. 



II) REVIEW OF THE EXACT SOLUTION OF 
THE BCS MODEL 



The red 
Hamiltonia: 



BCS model is defined by the 



H B cs= £ i° c )<y c 3° - 9 d J2 c ]+ c ]- c J'~ c 3'+ W 

3,<r=± 3,3' 

where Cj t ± ( resp. c ] ±) is an electron destruction ( 
resp. creation) operator in the time-reversed states \j, ±) 
with energies £j, d is the mean level spacing and g is 
the BCS dimensionless coupling constant. The sums in 
(|l|) run over a set of f2 doubly degenerate energy lev- 
els Ej(j — 1,...,0). We shall assume in this paper 
that the energy levels are all distinct, i.e. e< ^ Ej for 
i =/= j. The Hamiltonian (|l]) is really a simplified ver- 
sion of the reduced BCS Hamiltonian where all couplings 
have been set equal to a single one, namely g. This is 
the model that is commoly used to described ultrasmall 
grains, which describes the scattering of pairs of electrons 
between discrete energy levels that come in time-reversed 
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states. Hereafter we shall refer to ([!]) simply as the BCS 
Hamiltonian. 

As mentioned in the Introduction, Richardson had long 
ago solved this model exactly fori-aja-arbitrary set of lev- 
els, Sj, not necessarily all distinctJiiTE-!. To simplify mat- 
ters, we shall assume that there are not singly occupied 
electronic levels. As can be seen from (Q), these levels 
decouple from the rest of the system; they are said to be 
blocked, contributing only with their energy Sj to the to- 
tal energy £. The above simplification implies that every 
energy level j is either empty ( i.e. |vac)), or occupied by 
a pair of electrons ( i.e. ct + ct _|vac)). Denote the total 
number of electrons pairs by N. Then of course N < f2. 
The most studied case in the literature corresponds to 
the half-filled situation, where the number of .electrons, 
N e = 2N , is equal to the number of levels fflaPJl In the 
absence of interaction ( i.e. g — 0), all the pairs occupy 
the lowest energy levels forming a Fermi sea. The pair- 
ing interaction promotes the pairs to higher energies and 
eventually, for large values of N, all the levels are pair 
correlated, giving rise to superconductivityfil. 



Richardson's solution 

In order to describe Richardson's solution one defines 
the hard-core boson operators 
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which satisfy the commutation relations, 



The Hamiltonian ([j]) can then be written as 



(2) 
(3) 
(4) 



where we have set d = 1 ( i.e. all the energies are mea- 
sured in units of d). Richardson showed that the eigen- 
states of this Hamiltoniapqwith N pairs have the (un- 
normalized) product formEJltLS ( for a [direct proof of the 
results of this subsection see referencaHj ) 



l^>fl = IliLi^|vac), B v 



(5) 



where the parameters e v (u = 1, . . . , N) are, in general, 
complex solutions of the N coupled algebraic equations 



The ground state of Hbcs is given by the solution of 
eqs.(|]) which gives the lowest value of £ (JY). The (nor- 
malized) states (||) can also be written asE^ 

C 



\N)i 



i> R {h,---,3N)bl---b] N \™) (8) 



where the sum excludes double occupancy of pair states 
and the wave function ip takes the form 



A' 



•^)=eii 



i 



V k=l n 



(9) 



The sum in (|9|) runs over all the permutations, V, of 
1, • • • , N. The constant C in (|J) guarantees the normal- 
ization of the stateE3 ( i.e. r(N\N)r = 1); its expression 
will be given in section IV. 

A well known fact about the BCS Hamiltonian is that 
it is equivalent to that of a XY model with long range 
couplings and a "position dependent" magnetic field pro- 
portional to Sj. To see this let us represent the hard-core 
boson operators (g) in terms of the Pauli matrices as fol- 
lows, 



in which case the Hamiltonian (jij) becomes 



H 



BCS 



H 



) ~ j "i ~ 



H XY = -j: i 2e j q-i(T+T 



g(0/2 -iV) 



where the matrices 



T a = T.%it a j (o = 0,+,-) 



(10) 



(11) 



(12) 



t° = ±<7 Z 



tj = cr ■ 



satisfy the SU(2) algebra, 



[T a ,T b ] = f ( c lb T b 

ft = f°r = -i, Jo + ~ = 2 



(13) 



whose Casimir is given by 



-(T + T~ 
2 V 



+ T-T+) 



(14) 



N 

E 



E 



i 



25, 



(6) 



The energy of these states is given by the sum of the 
auxiliary parameters e„, i.e. 



Integrability of the BCS Hamiltonian 

From the existence of an exact analytic solution of 
Hbcs, one may expect that Hbcs shouldr-be-integrable. 
Indeed CRS found the integrals of motionc§Ej, 



N 



(7) 



Ri = -%-9 E 



tj • tj 



, (* = i,...,n) (is) 
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where the denominator does not blow up since we are 
assuming non degenerate energy levels. Integrability 
amounts to the eqs. 

[H B cs,Ri} = [Ri,Rj] = 0, {i,j = (16) 

Denote the eigenvalue of Rj acting on the state (||) by 
Xj, namely 



R j \N) R = X j \N) R 



(17) 



CRS, seemingly unaware of Richardson's solution, did 
not give an expression of Xj in their work. However, they 
did show that Hxy given in eq.(|Tl|) can be expressed in 
terms of the operators Ri as 



H XY =Y d ^ j R j +gC£ d R j f 



(18) 



Hence, combining (0) and (HI) one can find the eigen- 
values of Hxy 



Sxy = Y, 2£ ^ + 9{Y, X j) 2 -- a 9^ 



(19) 



and in turn those of £bcs by recourse to eq.(|Tl]). We 
shall show in section IV that Xj has the simple expression 




a 

- y 

J=l(*i) 



(20) 



One can check this result by deriving the energy (0) from 
eqs. (0), © and ©. 

This ends the presentation of the exact solution of 
the BCS Hamiltonian. The existence of an underly- 
ing analytic structure reminiscent to that of a CFT is 
apparent from eqs. (||), ( |l5| ) and (|2(]). Indeed, the 
aforementioned equations contain factors of the form 
l/(z — z') where z and z' stands for either 2Ej or e„. 
Terms of this sort arise quite naturally as correlators ( 
i.e. (A(z)B(z')) = l/(z—z')) of chiral primary fields A(z) 
and B(z) in diverse CFT's. The problem is to identify 
which CFT explains all the features presented so far in a 
unified manner. We shall argue that the solution of this 
problem is given by the Wess-Zumino-Witten (WZW) 
model based on the affine Kac-Moody group SU{2)k in 
the limit where the level k goes to —2. The proof of this 
result requires standard tools of CFT and, more precisely, 
the free field or Coulomb Gas (CG) representation of the 
WZW model. 



Ill) REVIEW OF THE FREE FIELD 
REPRESENTATION OF THE Sf/(2) A -WZW 
MODEL 



The material presented in this section is standard in 
the CFT literature^. Nevertheless, we have included it 



for the benefit of readers that are not experts in CFT. 
This will allows us to highlight the main tools we shall 
use in later sections. We shall follow closely referenced . 

The WZW model is an interacting theory which nev- 
ertheless admits a description in terms of free fieldsE3. 
The correlators and conformal blocks can then be easily 
calculated as integrals of vacuum expectation values of 
vertex operators. This gives an integral representation 
of the conformal blocks which satisfy automatically the 
Knizhnik-Zamolodchikov (KZ) equationsEj. In the case 
of the SU(2) k -WZW model the free fields are a j3 - 7 
system with conformal weights 1 and 0, respectively, and 
a boson field ip which |Saiiafy the following operator prod- 
uct expansion (OPE)UH, 



tf(z) ip(w) — — ln(z 



(21) 



w 



The WZW currents J a (z)(a — 0,±) can be expressed 
in terms of these fields as ( hereafter, normal order of 
operators will be implicitely assumed) 



J+ = i /3 
J = i [h 2 + ±ld^~ fcd 7 ] 



(22) 



which satisfy the OPE's 



J a (z) J b (w) 



k/2 
(z—w) 2 



1 fab rc 
z-iii Jc J 



reg. terms (23) 



where f® b are the SU(2) structure constants defined in 
eq.@, and q 0Q = 1, q+~ = q~+ = 2. The level of the 
WZW model, k, is related to the "charge" ao by the eq. 



k + 2 = 



1 



(24) 



If k is a positive integer, the WZW model is a Ratio- 
nal Conformal Field Theory (RCFT) with k + 1 primary 
fields labelled by the total spin, j = 0, 1/2, . . . , k/2. In 
our CFT approach to BCS, we shall need to consider 
the limit where k — * —2, which corresponds to taking 
ao — > 00. This is a singular limit which takes us away 
from the rational WZW models. Actually, the case when 
k is exactly —2 is mathematically interesting due to its 
relation to the singular hyperplanes in the representa- 
tion theory of affine Kac-Moody algebrasaH. For non- 
positive integer values of k, we can still define the theory 
by the free field representation given above. 

The Sugawara energy-momentum tensor T$ U g of the 
WZW model is given by the sum of the energy momen- 
tum tensors of the (3 — 7 system and the bosonic field 



T S ug ^ f3d 1 -\(d V ) 2 +ia 8 2 <p 



(25) 



The central extension c of the Virasoro algebra gener- 



ated by the modes of Ts ug (= J2 n L 



is 



3 



c = 3 - 12aS = 



3k 



(26) 



where the (3 — 7 system contributes with 2 and the field 
ip contributes with 1 — 12q;q. In ( |26| ) we have used the 
relation (^), which for integer k's gives the well known 
value of the Virasoro central charge of the SU (2)fc-WZW 
model. In the limit (k + 2) — * the central extension c 
diverges. In order to get a meaningful theory one has to 
scale the Virasoro operators as L n = limk^-2(k + 2)L n . 
In that limit the Virasoro algebra becomes 



[L n , L m ] 







(27) 



which suggests some sort of integrability. In fact the 
commutativity of the Virasoro operators L n has been 
used to study the representation theory of the SU(2)k=-2 
Kac-Moody-, algebra for the proof of the Kac-rKazhdan 
conjectures concerning character formulas (seet^l for ref- 
erences). 

The primary fields & m (z) of the WZW model are la- 
belled by the total spin j = 0,1/2,... and the third 
component of the spin m — j, . . . , —j. Their free field 
representation is given by 



(28) 



V a Az) = e iQ ^W, 



-2a j 



and have a conformal weight Aj given entirely by that of 
the vertex operator V a , namely 



1) 



(29) 



In the free field representation of a CFT, every pri- 
mary field has a conjugate version besides its "direct 
representapt." , which is needed for the computation of 
correlators^. In the case of the WZW model, the conju- 
gate of the primary field with m = j is 

^(z) = (3 s +^(z) V 2ao[s+3) {z), s = -(k + l) (30) 

The corresponding equation for m < j is much more 
complicated, and it is a sum of terms where the difference 
between (5 fields and 7 fields is given by s + j + m. 

A particular case of ( pfj| ) is when j = 0, which corre- 
sponds to the conjugate field of the identity 



I(z)=f3 s (z) V 2aoS (z) 



(31) 



A consequence of (^Tj) is that the expectation values of 
operators should satisfy the following charge neutrality 
conditions 

Np-N 7 = s, ^a 4 = 2a s, s = -(k+l) (32) 



where Np and N 7 is the number of (3 and 7 fields in the 
correlator, and ot^ are the charges of the vertex operators 
made of the field <p. Eq.(|32|) means that there is a back- 
ground charge — 2aos in the boson sector and a charge 



— s in the (3 — 7 sector, which need to be neutralize for 
the correlator to be non-vanishing. The latter properties 
can alternatively be attributed to the out vacuum which 
have charges — 2ao and —s in the ip and f3 — 7 sectors, 
respectively. 

The remaining ingredient of the free field representa- 
tion is provided by the so called screening charge 



Q 



duS{u), S{u) = P(u) V 2ao (u) 



(33) 



whose basic property is that it commutes with the SU(2) 
current algebra and the Virasoro operators. In eq.(|33J) 
and below, du is meant to contain the factor l/(27ri) 
to takes care of the factor 27ri that comes out in the 
residue formula. Using the vertex representations (SM) 
and (|3^) of the primary fields, together with the screen- 
ing charge ([33]), one can compute the conformal blocks of 
the WZW model. Conformal blocks are the chiral build- 
ing blocks of correlators. The latter are obtained by com- 
bining the holomorphic and the anti-holomorphic confor- 
mal blocks and imposing monodromy invariance. In the 
WZW model a conformal block ip wzw (z±, . . . , Zn+i) in- 
volving f2 + 1 primary fields {(j'fc, inserted at 
the positions {zk}^l, can be associated with the SU(2) 
tensor product decomposition 



3i 



En 
k=i m k 



® jn -> jn+i 
mn+i 



(34) 



where (in+i, mu+x) appears as an outgoing state. The 
free field expression of the conformal blocks is given by 



(^ Zl ) . . .^-Jzn) ^Uzn+i) 
x § c du\S{ui) . . . § Cn duNS{u N )) 



(35) 



where <&mn+i is the conjugate of the outgoing state 
(jn+i, mn+i), and the screening charges are integrated 
along the contours C\, . . . , Cjy. The charge neutrality 
conditions (B2h applied to (B5|) yield 



jn+i, mu+i 



fc=i 



a 

E 

fc=i 



TOfc 



(36) 



which agree with the Clebsch-Gordan decomposition 
(|34|). The case when V = corresponds to the maximal 
allowed value of jn+i = S/b=iJfc- O n the other hand, 
if fl is even, the minimal value of ja+i is zero, which 
requires V = X)&=i 3k screening charges. Hence, the dif- 
ferent choices of the screening charges and contours give 
rise to all possible conformal blocks. In this manner the 
free field representation provides integral solutions of the 
KZ equations satisfied by the conformal blocks ( |35| ) . The 
KZ eqs areE3 



4 



Q+l 



K dz, ^ 



ip {zi, . . . ,2n+iJ (37) 



where k = (fc + 2)/2 and tj are the SU(2) matrices in the 
ji representation acting at the i th site. 



IV) CFT REPRESENTATION OF THE EXACT 
SOLUTION OF BCS 

Our first aim is to obtain Richardson's wave function 
(^) using the free field representation of the WZW model. 



Richardson's wave function 

The starting point is the speudospin version of the BCS 
model introduced in section II, according to which an 
empty energy level Si{i = 1, . . . , fi) has spin m, = 1/2 
while and occupied level has spin m, = —1/2. This sug- 
gests to rewrite Richardson's wave function as 



N 



4>, 



l (z 1 ,...,zn;e 1 ,...,e N ) = ^ JJ 



1 



V fe=l 



(38) 



where = 2si and Ik is defined by the condition mi k — 
— 1/2. The SU(2) quantum numbers m< of (p8|) satisfy 



E 1 

i=l 



= — - N 
2 



(39) 



Let us compare now the conformal block (|35|) and the 
wave function (|3q ). If we take jk = 1/2 (fc = 1, . . . , f2) in 
(Bq) and use (|39[), we are lead to the identifications 



n 



(40) 



which requires TV < tt/2. Hence from a formal point 
of view, we can regard Richardson's wave function as 

1 /2 

a conformal block involving ft primary fields ^ 
located at the positions Zj = 2ej, and a primary field 
^mnli ( z n+i) whose position we shall place at oo. The 
last ingredient we need, in order to reproduce Richard- 
son's wave function using CFT tools, is to find the role 
played by the BCS coupling constant g. We shall see 
below that g is associated to the operator 



V g = exp(- 



9 



dz zdip(z)) 



(41) 



Our claim is that Richardson's wave function (|3^) is 
given, up to a proportionality factor, by the limit 



^(z,e)cx lim Vm G (z) 

OtQ— >oc 



(42) 



where the Coulomb Gas wave function ip CG is given by 
the following expectation value 



= {Vg^izx) . . .®?n n {zn) 
x *&(°°) § Cl duiS( Ul ) ...§ Cn du N S{u N )) 



(43) 



Except for the presence of the operator V g , eq. ( f43| ) coin- 
cides with the conformal block (|5|). Let us now prove 
cq. ji^ ) . First of all, using the free field representation of 
the primary fields, one can write (^) as 

C G (z) = § Cl du i--- §c N duN^iz, u) ^(z, u) (44) 
where 

w = <^n-Li^ Q o(^)n!Li^oK)^ QO ( S+Jn+1 )(oo)) (45) 



N 



Ip. 



01 _ 



Using eqs.(^l|) and the Wick theorem, one can show that 
ipP J is, up to a sign, equal to the Richardson's wave func- 
tion, namely 



V£(z,e) = (-l) N ^(z,e) 



(46) 



In this equation we have choosen the positions of the 
screening operators Ui e qua l to the Richardson's param- 
eters ej. However, in eq.(|44|) one must integrate over the 
it's. The job of the limit ctQ — > oo is to set Ui — e^. Let 
us see how this happens. First of all the contribution of 
the vertex operators can be computed using the formula 



{V g V ai (w 1 )...V aM (w M )) 



(47) 



where the contour C g encircles all the coordinates Wj. 
The factor tp v then becomes 

■0 (Z, U) = J] e a o z ^9 e -2ag«v/9 



(48) 



IW 



r° n,<>* 



Vinr, 



n, 



-2<.>; 



In the limit ao — > oo, the integral (|4J) can be computed 
using the saddle point method. Indeed, writting ([48]) as 



^(z,u) 



-ag C7(z,u) 



(49) 



where 



U = - Ef<i m fe - z j) - 4 E^<u m K 



g 1 



(50) 



the stationary solutions of U are given by the solutions 
of Richardson's eq.ra), namely 



5 



3u„ 



1 N 

-+ E 



E 7~7~ ' 



7/ — p 5 i / / \ °t / 1 ^7 °^ 



(51) 

Under these conditions the saddle point value of the in- 
tegral (El) is 



-agff(.,e) ^7( Z)G ) (52) 



(27T 1 /2 Qo )iV( dctA )l 

where A is the N x N hessian matrix defined as 



au 



A, 



E 



i (z;— e„) 2 



e: 



2 



(53) 



A 



Eqs.©, @ and (f|) constitute the proof of @. 
Moreover the factor (detA) -1 / 2 in ( p32] ) turns out to co- 
incide with the normalisation constant C appearing in 
the normalized state (||)c3, namely 



C= l/(detA) 



1/2 



(54) 



coupling version of BCS, which is valid for all val- 
ues of the BCS coupling constant g, since there are 
not phase transitions as we go from weak to strong 
coupling. 

• So far, we have assumed for simplicity that the 
pair degeneracy of the energy levels £j is unity, 
which corresponds in the WZW model to primary 
fields with spin jk = 1/2. The construction can be 
straighforwardly generalized to levels with higher 
degeneracy, dk — 1,2,..., in which case the asso- 
ciated primary fields have spin jk t~ \dtJ2. I n this 
case the Richardson equations readE^Ea, 



1 



N 

E 



E 



(55) 



and they can be derived using the free field rep- 
resentation explained above in terms of the fields 



$# fc (z fc ) with j fc =4/2. 



Let us make some comments about the results pre- 
sented so far. 

• Richardson has observed that eqs.(||) could be de- 
rived as the stationary configurations for a set of 
electrostatic charges in 2Do. The electrostatic po- 
tential of these charges is given by a^U or rather 
by the sum a 2 , (£7 + U*). Our CFT derivation of 
the exact BCS solution shows that the value of the 
charges are — «o for the energy levels, 2ao for the 
screening ones, and a charge 2aa(s + jn+i) placed 
at infinity. The sum of all these charges neutral- 
izes the background charge — 2aoS. The stationary 
conditions are set by the limit ao — > oo, which is 
equivalent to the limit fc — > — 2. 

• The operator V g , which is needed in order to get 
the term 1/g in the Richardson's equations, breaks 
conformal invariance in an explicit manner. This 
is why ip CG is not, strictly speaking, a conformal 
block of the WZW model—In the spirit of Perturbed 
Conformal Field TheoryED we could say that the 
WZW model has been perturbed by the chiral op- 
erator — ^ <f> dzzdip which is equal to — a a_i/g, 
where a_i is the n = — 1 mode of the field dtp 

(i^W = E„^- n - x ). 

• The perturbative renormalization group (RG) anal- 
ysis of the BCS model yields that the coupling con- 
stant g > flows to large values in the infrared 
regime, leading lo the superconducting instability 
of the Fermi seaLj. In our approach the fixed point 
of this RG flow is described by the SU{2)k=~2- 
WZW model. The 1/g perturbation takes us away 
from this fixed point, breaking conformal invari- 
ance. In that sense, we are dealing with the strong 



The integrability of BCS from CFT 



At this stage it is clear that eqs. (|15|) and ( |17| ) must 
be related to the KZ eq. (p7|). After all, in the limit 
ao — > oo the CG wave function (|4|) coincides with the 
WZW wave function. Let us show this in detail. The 
KZ eq. (|37|) with zn+i — oo and the definition Jl5| ) can 
be written as 



Let us define a new wave function if) 







if, 



WZW _ H xy /2 9 k 



(56) 



(57) 



Using the commutativity of Ri and Hxy and the fact 
that 



d_ 



Hxy — —t { j 



eq.(p9) becomes 



1 d . 

— Ri + k— V 

2g OZi 







(58) 



(59) 



This equation is completely equivalent to the KZ eq. 
Imposing that if) diagonalizes Ri with eigenvalue Xi , if) be- 
comes proportional to the Richardson solution, or rather 
to if> CG . Hence eq.([l9|) yields the following expression for 
\i 



1 ~ 2 d Zl 



(60) 



which in turn leads to the formula (|20|). In other words, 
if> wzw , ip CG and ip R agree up to overall factors in the 
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limit when ao — > oo. It is interesting to observe that the 
factor relating the WZW and Richardson wave functions 
has indeed a conformal structure reflected solely in the 
presence of logarithms, namely 

^WZW = e *l(lE XY -U) ^fl ( 61 ) 

The Coulomb energy U contains in a disguised manner 
both the Richardson's eqs. and the constants of motion 

V) THE BCS/CFT CORRESPONDENCE 

In the previous section we have shown the closed rela- 
tionship between the BCS and the CFT wave functions. 
In this section we want to investigate this relation at 
the second quantized level. In a Quantum Field Theory 
a generic N-body wave function ^(xi, ■ ■ ■ , xn) is usu- 
ally constructed from the overlap of a N— body state \$f) 
with the eigenstates \x\, . . . , xn) created by the action of 
the field operator ip(x) acting on the Fock vacuum |0), 
namely 

*(asi, . . . , ajjy) = <*|an, . . . , arjv) (62) 
|aii, — , aciv) =i){x\)...'ip(x N )\Q) 

An interesting example of this formalism is provided by, 
the CFT interpretation of the Laughlin wave functionEZl 
of the Fractional Quantum Hall effect (FQHE), first pro- 
posed by Fubinill. In that case tJj{x) is a vertex operator 
of a single boson governed by a c = 1 CFT. One can also 
make use of an array of bosons. 

In the spirit of the FQHE/CFT correspondence, one 
can interpret the Coulomb Gas wave function (|4^) as 

V'm ( z ) = \zi,mi, . . .,za,ma) (63) 

where 

|zi,mi, . . .,zn,ma) = nf=i^ i (^)|0) (64) 
(^ GG \ = (oo(n - 2N)\V g n£Li § 0u du„S(u u ) 

The out vacuum (ao(fl — 2N)\ has been defined as the 
action of the operator $mn+i on the out vacuum of the 
WZW model, and has a charge ao(£l— 2N) in the ip sector 
and no charge in the (3 — 7 sector. The different states 
(ip CG \ correspond to different choices of the integration 
contours C v of the screening operators S(u„), which in 
CFT yield different conformal blocks. For example, the 
ground state is obtained by choosing the z/ h -contour {y = 
1, . . . , N) to run from z v to infinity, where we assume 
that the energy levels are ordered in increasing order, 
i.e. z\ < z 2 < ■ ■ ■ < Zq. The excited states correspond to 
other contour choices. The total number of eigenstates of 



Hbcs with N pairs and £1 levels, which is given by the 

combinatorial number f ^ J , coincides with the total 

number of contour choices. 

Eqs.(|63|) and (|64|) provide the basic correspondences 
between the exact solution of BCS and CFT, which can 
be extended to other instances. We collect them in table 
1 and make some comments below. 



BCS 


CFT 


Pair energy level 


WZW Primary field 


Pair degeneracy (dk) 


Total spin (j k = d*/2) 


Eigenstates of Hbcs 


Conformal blocks 


Richardson's eqs 


Saddle point conditions 


Integrability (CRS) 


KZ eqs. 


g = 00 { g finite) 


WZW ( Perturbed WZW) 


Cooper pair operator 


Screening operator 


Phase stiffness 


cx ag 


empty, occupied level 


spin up, down 


nf=i & li vac > 


nlLi< 2 te)io) 


r(n\ 










§ Zj dzJ+(z),J-(z),J°(z) 


C 


1/VdetA 


A, 


9 dU c 

2 dzt 






Hxy 


-9^0 



Table 1. The BCS/CFT correspondence 



• In analogy to eqs.(|6J) one could try to define the 

bra state(zi, mi, . . . , zn, mn\ using the conjugate 
~ a 

operators $, 2 n (z). However, this state would have 
a large charge which would be difficult to compen- 
sate for. Similarly, the ket state \ip GG ) is not easy 
to construct since the screening operator S(z) does 
not have a conjugate version. A possibility would 
be to use the second screening operator, given by 
S(z) = /? s-1 exp(— iip(z) / a p3. This operator ap- 
pears in the construction of fractional spin repre- 
sentations which are related upon a certain reduc- 
tion to the minimal models. In the limit ao — > 00 
we see that S(z) converges to the identity. Further 
work is needed to clarify this issue. 

• In section III we defined the Virasoro operators L n 
as the singular limit limfc^_2(fc + 2)L n . From eq. 
( p9| ) we see that Ri could be identified with the ac- 
tion of — g_L-\ on ip GG , where L_ x = d/dzi. Simi- 
larly, eq.(|f8|) implies that Hxy could be identified, 
up to constants, with —g(k + 2) J2j z j9/dzj and 

thus with —qLq. Upon these identifications, the 
integrability of the BCS model, given by eqs . (|lq ), 
becomes equivalent to the commutativity (Mi of 
the Virasoro operators Lq and L_ t . We may expect 
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the existence of another CRS-like integrals of mo- 
tion associated to the Virasoro operators L n (n < 
-2). 

• The Richardson's state r(N\ corresponds, in the 
CFT formulation, to the state (ip CG \, which is the 
product of N screening operators acting on the out 
vacuum (ao(fl — 2N)\. This correspondence sug- 
gests that in the grand canonical (g.c.) ensemble ( 
where the number of pairs is not fixed ) the corre- 
sponding state should be given by 

(V£ c G | = (0|ex P (^dzS(z)) (65) 

( Note that we have assumed the half filled condi- 
tion tt = 2N1 This state is similar to the BCS 
state given bytl 

{BCS\ = (vac| exp(^ 5fe 6 fc ) (66) 

k 

where g k is the ratio Vk/uk of the BCS variational 
parameters. In this sense, the screening charge 
§ dzS(z) is the CFT version of the Cooper paki 
operator ^2 k gkbk- In CFT it has been arguedES 
that the screening operators can be exponentiated 
into the action and that their number is fixed upon 
imposing the charge neutrality conditions on cor- 
relators. This CFT exponentiation corresponds to 
working in the grand canonical emscmble in BCS. 

• In the previous item we argued that S(z) — 
P(z)exp(2ao'np(z)) is the CFT analogue oigtbk- On 
the other hand, bk corresponds to the contour in- 
tegration of J + (z) — i/3(z) around the pair energy 
Zfe = 2efc. Hence it is natural to associate the BCS 
variational parameter g/. with exp(2aoi<p(z)). This 
means that 2ao(p(z) can be associated to the phase 
of the superconducting order parameter. Shifting 
ip(z) by a constant leads to an overall phase shift 
of the BCS order parameter. This correspondence 
yields an insight about the physical meaning of q.q, 
which seems to be related to the phase stiffness or 
the superfluid density n s . Indeed, if we identify 
the phase of the superconducting order parameter 
9 with ao<p then the Lagrangian of tp becomes that 
of a continuum XY model for 9, with Oq playing the 
role of the superfluid density. Actually appears 
in the denominator of the Lagrangian while n s ap- 
pears in the numerator. However recall that we are 
working in the energy space so things are inverted. 
The identification of 9 with aoip is also consistent 
with the fact that both variables are defined mod- 
ulo 2n. The limit ao — > oo therefore corresponds to 
the limit of very large phase stiffness which in fact 
leads to the standard BCS theory, where the phase 
of the superconducting order parameter is rigid and 



plays no role in fixing the critical temperature pr. 
other observablesEj. As was shown by RichardsonEfl 
cqs.^) reduce in the bulk limit N — > oo to the BCS 
gap equation and hence the state (0) becomes the 
fixed N projection of the mean field g.c. BCS state. 
Finite values of ao should lead to non mean field 
theories with the phase 9 playing a dynamical role. 
It is rather intriguing that models of this sort have 
already been proposed by several authors for an 
explanation of high-T c superconductivityo. 



VI) COMPARISON BETWEEN BCS AND THE 
FQHE 

In sections V and VI we noticed some analogies be- 
tween the CFT approaches to BCS and the FQHE. Let 
us consider them in some more detail. A common feature 
is the Coulomb Gas treatment. In the FQHE the CG is 
associated to the Laughlin-wave function of N e electrons 
at filling factor v = l/m,E3 

■ ■ ■ ,w Ne ) = Y[(wi - Wj ) m e-^ l^l 2 (67) 

i<j 

The norm of (pTj) can be seen as a classical probability 
distribution e~^^ of a two-dimensional one-component 
plasma at ficticious temperature [3 — X/m and potential 
energy Ul wherdf3, 

N e 

U L = -2m 2 Y J Mw l -w 3 \ + 7 ^Y.\ w ^ ( 68 ) 

The particles with charge m repell each other with a log- 
arithmic interaction, and they are attracted to the ori- 
gin by an uniform neutralizing background charge with 
density p — 1/(2tt£%), where £b is the magnetic length, 
which has been set equal to one in ( |67| ) and (p8|) . For 
small values of m(= 3, 5, . . .) the electrons form a liquid 
with uniform density p e = l/{2 r r:ml 2 B ) which neutralizes 
the background charge. However, for large values of m, 
Quantum Montecarlo studies have shown that the Laugh- 
lin liquid becomes a solid ( i.e. a Wigner crystal) where 
the positions of the charges are localized. 

The comparison between the wave function ip v given 
in eq. ( ff8| ) and the Laughlin wave function ipL suggests a 
formal identification of the electron positions uij with the 
screening positions Uj rather than with the pair energies 
levels Zj. The reason is that both the u's and the w's 
are subject to integration, while the z's are held fixed. 
Following this analogy, we may stablish the relations 

k + 2 

m = 4q;q, or v = — - — (69) 

according to which the freezing of the screening charges 
in the limit ao — > oo would parallel the Wigner crytal 



8 



structure of the FQHE when m is large. On the con- 
trary, for finite values of ao the screening charges, which 
are essentially Cooper pairs, would delocalize becoming a 
sort of liquid. The discussion at the end of the previous 
section suggests that this liquid should arise from the 
fluctuations of the phase of the superconducting order 
parameter. 

Besides these analogies, the Laughlin and BCS 
Coulomb Gas models differ in the nature of the back- 
ground charge. In the Laughlin case the charge is two- 
dimensional, while in the BCS case the linear terms ap- 
pearing in U ( see eq. (|50|)) can be attributed to a linear 
uniform density pi oc ao/g placed at infinity. The lat- 
ter density creates a uniform electric field d x (U + U*) 
along the x = (z + z*)/2 axis. Another difference is that 
the BCS theory is not really conformal invariant for fi- 
nite values of g while the Lausklin state has gapless edge 
excitations described by CFTc2l. 

VII) CONCLUSIONS AND PROSPECTS 



modelsBEl. For the G^-WZW model the singular limit 
is given by, k + h — > 0, where h is the dual Coxeter num- 
ber of (Jk3. The charge and spin independent pairing 
Hamiltonians studied by RichardsonE3Ej in the context 
on Nuclear Physics probably belong to this category of 
models. 
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Note added After completion of this work we have 
been informed by prof. JL Belavin about some related 
work by H.M. Babujiano, where he applies the Bethe 
ansatz and the Jinizhnik-Zamolodchikov eqs. to the 
Gaudin magnetsEl The Gaudin's model is given essen- 
tially by the g — > oo limit of the reduced BCS model. 
In fact the Gaudin's Hamiltonians can be identified with 
lim^oo Rj/g, where Rj are the CRS conserved quanti- 
ties define in eq. (|l5|) . 



In this paper we have established a closed relationship 
between the exact solution of the BCS Hamiltonian and 
the Coulomb Gas version of the SU(2)k-WZW model in 
the singular limit when k — > —2. The Richardson's wave 
function comes from the (3 — 7 chiral correlators, while 
the Richardson's eqs. and the normalization factor of 
the state arises from the saddle point evaluation of the 
chiral boson correlators. The BCS coupling constant g 
enters the construction as a perturbation of the WZW 
model and breaks conformal invariance. The integrabil- 
ity of the BCS model is related to that of the WZW 
model through the KZ equations, which has lead us to 
an expression of the integrals of motion of the BCS model 
found by Cambiaggio, Rivas and Saraceno. We have pro- 
posed a BCS/CFT correspondence which, in many re- 
spects, parallels the CFT interpretation of the Fractional 
Quantum Hall effect. We have conjectured that the sin- 
gular limit a§ — > 00 amounts in physical terms to the 
limit of very strong phase stiffness, which leads to the 
mean field BCS theory. Finite ao generalizations of the 
BCS model may correspond to non mean field theories 
where the phase of the superconducting order parame- 
ter should have a dynamical role as in some models of 
high-T c superconductivity!^]. 

Besides giving new insights into the exact solution of 
the BCS model, the CFT approach may also help in 
solving some problems as the computation of observables 
with the Richardson's exact solution. The finite temper- 
ature BCS model is also an interesting problem which 
one may try to address with CFT tools. 

Finally, the BCS/CFT approach can be straighfor- 
wardly generalized to any WZW model based on an 
affine Kac-Moody algebra Gk, where k is the level and 
G is a semi-simple Lie group or supergroup. As in the 
SU(2) case, one can use the free field realization of these 
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